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PACS 78 . 67 . -n - Optical properties of low-dimensional, mesoscopic, and nanoscale materials and 
structures 

PACS 73.20.Mf - Collective excitations (including excitons, polarons, plasmons and other charge 
density excitations) 

PACS 71.45.Gm - Exchange, correlation, dielectric and magnetic response functions, plasmons 

Abstract. - We propose an explanation for the anomalous behaviour observed in the relaxation 
dynamics of the differential optical transmission of noble-metal nanoparticles. Using the temper- 
ature dependences of the position and the width of the surface plasmon resonance, we obtain a 
strong frequency dependence in the time evolution of the transmission close to the resonance. In 
particular, our approach accounts for the slowdown found below the plasmon frequency. This in- 
terpretation is independent of the presence of a nearby interband transition which has been invoked 
previously. We therefore argue that the anomaly should also appear for alkaline nanoparticles. 



Absorption and transmission of a laser beam by a metal- 
lic nanoparticle unveil the properties of its conduction 
electrons. The use of femtosecond lasers in pump-probe 
spectroscopy then gives access to the electron dynamics 
on extremely short time scales [1-6]. In particular, the 
time dependence of the transmission spectrum during and 
after the excitation by the pump laser allows to follow the 
details of the electron relaxation. 

In pump- probe experiments, a pump laser with a wave- 
length much larger than the diameter of the nanoparticle 
couples to the centre of mass of the electron gas, thereby 
exciting the surface plasmon mode. The collective ex- 
citation decays on a very short timescale of the order 
of 10 fs. Subsequently, electron-electron interactions lead 
quite rapidly to the thermalisation of the electronic system 
at an elevated temperature. Only on longer timcscalcs, of 
the order of 100 fs up to picoseconds, the electron-phonon 
coupling leads to an equilibration of the electronic system 
with the lattice degrees of freedom. In the case of noble- 
metal nanoparticles, this picture has to be completed by 
the possibility of interband transitions which are absent 
in alkaline nanoparticles. 

In ref. [2] the relaxation of the differential transmission 
spectrum of copper nanoparticles was found to exhibit a 



slowdown close to the plasmon frequency. This feature 
persists far beyond the lifetime of the plasmon excitation 
and therefore was attributed to the slowdown of the en- 
ergy transfer from the electrons to the lattice for frequen- 
cies close to the plasmon resonance. These considerations 
lead to the following puzzle: How can the plasmon reso- 
nance frequency play a role for relaxation processes which 
occur on timescales much larger than the lifetime of this 
excitation? 

The s-d multiband transition, which for copper has its 
onset close to the plasmon resonance, has been invoked 
as a possible explanation of the unexpected slowdown. A 
many-body effect based on the interband resonance scat- 
tering of d holes with surface plasmons was proposed to 
lead to a strong frequency dependence of the relaxation [7] . 

Subsequent experiments with silver nanoparticles also 
found the slowdown in the transmission dynamics [4]. 
Firstly, these experiments demonstrated the generality of 
the effect. Secondly, they established that the slowdown 
is only found above an excitation threshold. Finally, be- 
cause the interband transition in silver is far away from 
the plasmon frequency [6] , the explanation put forward in 
refs. [7] does not apply. 

In the present paper, we propose a generic explanation 
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of the observed frequency dependence of the transmission 
dynamics which does not rely on interband transitions. 
Therefore it applies not only to copper, but also to silver 
and even alkaline nanoparticles. Our mechanism is based 
on recently derived results for the finite-temperature cor- 
rections to the width and the position of the plasmon ab- 
sorption peak [8]. We show that the temperature depen- 
dencies of the two quantities lead to an anomaly in the 
evolution of the differential transmission spectrum with 
decreasing temperature. In particular, a slowdown in the 
relaxation of the optical transmission appears close to the 
plasmon resonance. 

The quantity that is measured experimentally is the 
frequency-resolved time-dependent transmission T on after 
the excitation of the system by a pump laser pulse. Then, 
the static equilibrium transmission %g of the sample mea- 
sured in the absence of the pump is subtracted to obtain 
the relative differential transmission [2,6] 
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Experimentally, the observed exponential decay of AT/T 
with time allowed to extract a relaxation time which 
as a function of frequency presents the above-mentioned 
anomaly. 

For the theoretical discussion we assume that after an 
initial transient time of about 100 fs the electronic system 
is thermalised with a time-dependent temperature T(t). 
In thermal equilibrium, the differential transmission can 
be related to the temperature-dependent absorption cross 
section a(u),T) by [9] 
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where T is the temperature of the nanoparticle before 
the excitation and a is its radius. This relation between 
transmission and absorption holds when the reflectivity of 
the sample is negligibly small, as it is the case in typical 
experimental situations with a low density of nanoparticles 
embedded in a glassy matrix. 

Assuming the Breit-Wigner form for the absorption 
cross section a of a nanoparticle in the frequency regime 
close to the plasmon resonance, we write 



er(u>,T) = s(a) 
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Here, oj sp (T) is the temperature-dependent resonance fre- 
quency of the surface plasmon excitation, 7(T) is the 
temperature-dependent width of the resonance, and s(a) 
is a normalisation prefactor which depends on the radius 
of the nanoparticle. 

The knowledge of the temperature T(t) as a function of 
time would give access to the time evolution of the differ- 
ential transmission. In principle, T(t) could be obtained 
for specific systems from a two-temperature model for the 
heat transfer from the electronic system to the lattice [10]. 



However, we choose not to follow this route since our ex- 
planation of the anomaly relies only on the fact that T(t) 
is monotonically decreasing. The qualitative and generic 
features do not depend on the system parameters entering 
in the quantitative description but on the temperature de- 
pendencies of uj sp and 7. For temperatures much smaller 
than the Fermi temperature Tp we can write 



Wsp (T)=<-a;g) (X) 
7 (T)= 7 (°)+ 7 ( 2 ) ( 
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(4) 
(5) 

The dependence of the positive quantities (Jsp , 7 < -°\ Wsp 
and 7 1 - 2 ) on the system parameters is specified below 
(eqs. (9)— (11)) , but it is not needed for the following qual- 
itative discussion. 

Inserting the Breit-Wigner form (3) of the resonance 
into the expression (2) for the differential transmission, 
one can see that a temperature-dependent shift of the 
resonance frequency lu sp leads to an anomaly in the re- 
laxation dynamics of the differential transmission. At a 
given temperature, AT/T results from the difference be- 
tween two absorption curves which are shifted in frequency 
with respect to each other and have different widths. The 
frequency dependence of AT/T is obtained from eqs. (2)- 
(5) and shown in fig. 1 for various temperatures. For a 
given temperature, there exist two frequencies 2 where 
the differential transmission vanishes. In the limit of small 
T/Tp these frequencies can be approximated by 




AT/T is negative for frequencies below oj\ or above Wjj, 
and positive between these two frequencies. As the fre- 
quency uj\ is very close to u>^ , our assumption of a Breit- 
Wigner lineshapc which is a good description close to 
the resonance does not represent an important restriction. 
However, this is not the case for u>2, which is farther away 
from uip\ For the parameters chosen in fig. 1, uj\ lies at 
the far right of the shown frequency interval. There, the 
detailed form of the absorption curve will be important. 
It might for example play a role whether one assumes the 
Breit-Wigner form (3) or the quasi-Lorentzian that has 
been used in ref. [6] to fit the plasmon resonance over a 
large frequency range. In addition, for the case of copper 
the effects occurring around wJj will be masked by the in- 
terband transition. Therefore we focus our discussion on 
the anomalous behaviour close to 

When the temperature of the electronic system de- 
creases after the initial excitation, the position of the 
absorption peak moves back towards its low-temperature 
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Fig. 1: The differential transmission for ef/Hljm = 1, £f/Vo = 
0.5, kpa = 15, and To = (corresponding to a sodium nanopar- 
ticle with diameter 3.3 nm) is shown as a function of the fre- 
quency u). The temperatures are T/Tf = 0.05, 0.035, and 0.02 
for the solid, dashed, and dotted line, respectively. The inset 
shows an enlargement of the region around u>l where AT/T 
changes sign. 



frequency During this relaxation process, the differen- 
tial transmission is reduced in amplitude. The shift of 
u\ towards higher frequencies and the reduction in ampli- 
tude have opposite effects on the differential transmission 
for frequencies below u>\. The contributions subtract and 
lead to a slowdown of the relaxation. As an example, one 
may consider the crossing of the solid and the dashed line 
in the inset of fig. 1. Here, the differential transmission 
at T/Tp = 0.035 (dashed line) has the same value as at 
the higher temperature 0.05 (solid line). In contrast, for 
frequencies slightly above wf, the two effects add up and 
the relaxation is accelerated. A very special dynamics can 
occur at frequencies which are crossed by uj\ during the 
temperature relaxation. In such a situation, the differen- 
tial transmission goes through zero and changes sign as 
the temperature decreases (see the inset in fig. 1). 

For the example of alkaline nanoparticles, where only s- 
conduction electrons are relevant and no interband tran- 
sition is present in the frequency regime of the plasmon 
resonance, we can show in detail how the mechanism pre- 
sented above leads to the anomaly in the transmission re- 
laxation. For this purpose, we use the temperature de- 
pendencies of the resonance frequency and of its width, 
which have recently been calculated [8] for temperatures 
well below Tv. 

In the approach of refs. [8] and [11], one decomposes 
the Hamiltonian H = H cm + H IC \ + H c describing the 
system of TV electrons in a spherical positive jellium into 
three parts [12]. The surface plasmon is represented by 
the centre-of-mass part H cm — hu)y[b^b, where b and 
are the usual ladder operators of the harmonic oscillator. 
Its frequency ujm = wm(1 — Aout/A) 1 / 2 is reduced with 
respect to the bare Mie frequency com = (47rn c e 2 /3TO c ) 1 / 2 
due to the so-called spill-out effect which accounts for the 



fact that a fraction N out /N of the electrons is found out- 
side the jellium sphere [13]. In the Mie frequency, e, m e , 
and n c denote the charge, mass, and bulk density of the 
conduction electrons, respectively. 

The relative coordinates constitute an environment 
leading to the decay of the plasmon. Within a mean-field 
approximation, one can write the corresponding part of 
the Hamiltonian as H Te \ = $^ a £ac£c a , where the oper- 
ators cjj populate one-particle states with energies e a in 
the effective potential V(r). These degrees of freedom are 
coupled to the plasmon by 



afi 
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where A = (fim c w^/2iV) 1/2 and 
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is a matrix element between two mean-field states. <d(x) 
denotes the Heaviside step function. 

Approximating the self-consistent mean-field potential 
by V(r) = VoO(r — a), using scmiclassical techniques (for 
details, see refs. [8,11]), and working up to second order 
in T/Tp, we obtain the expansions (4) and (5) for the 
position of the plasmon resonance and its width, where 
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Here, we have introduced the auxiliary functions [8, 11, 14, 
15] 
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^2 f 

g 2 (x) = — {y/x(x + l)(6x - 1) + In (y/x + y/x + l) 



low-temperature Sommerfeld expansion of the electronic 
energy one gets the temperature [11] 



-e(a;-l) y/x(x-l)(6x + l) 

+ In (v/i + \/x - 1) }. 



£f, vf, and &f are the Fermi energy, velocity, and wavevec- 
tor, respectively. 

The dependence on x of £o and £2 arises from the cou- 
pling of the plasmon to the relative coordinates while the 
dependence on y is due to the spill-out effect. Both ef- 
fects depend quadratically on the temperature. The pa- 
rameter j] is a cutoff of order unity that appears because 
the low-energy particle-hole excitations form the plasmon 
and thus do not contribute to its damping [8, 15]. For 
numerical purposes we take r] — 0.5. It is important to re- 
mark that the finite temperature corrections scale as l/o, 
and therefore become negligible for sufficiently large clus- 
ters. For a quantitatively correct description of very small 
nanoparticlcs (a < lnm), nonmonotonic size-dependent 
corrections have to be added [16]. In the case of noble- 
metal nanoparticles embedded in a glass matrix, the effect 
of the d electrons and the dielectric mismatch at the bor- 
der of the nanoparticle should be taken into account [17]. 
In the present context, the essential feature of the func- 
tions 92 and £2 is that they are positive. Therefore higher 
temperatures lead to a redshift of the resonance and to an 
increase of its width. 

Numerical calculations within the temperature- 
dependent time-dependent local density approximation 
predicted a redshift of the resonance up to a system- 
dependent temperature, followed by a blucshift at higher 
temperatures [18]. This nonmonotonic feature has not 
been observed experimentally to our knowledge and 
is not captured by our low-tcmpcrature expansion. 
Systematic studies of the frequency shift and plasmon 
width have been carried out in gold [3] and silver [6] 
nanoparticles yielding a decrease of both corrections 
with increasing pump-probe time delay, i.e., decreasing 
electron temperature. 

In fig. 1 we used the analytic results presented in eqs. (4) 
and (5) for the temperature dependence of the plasmon 
resonance to describe the evolution of the differential 
transmission with decreasing electronic temperature. The 
curves present the qualitative features of the experimen- 
tally observed differential transmission of refs. [2] and [6]. 
As discussed after (6), the relaxation of the tempera- 
ture leads to an anomaly at frequencies close to the plas- 
mon resonance. We now investigate the properties of this 
anomaly in detail. 

For different frequencies below (solid lines) and above 
(dashed lines) the frequency lo\ we show in fig. 2 the re- 
laxation dynamics normalized by the value at the maxi- 
mum temperature Tm — 0.05Tp reached after the excita- 
tion. An estimate for Tm can be obtained by assuming 
that the nanoparticle absorbs n photons from the pump 
laser and increases its energy by about nfko sp . Using the 
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For the sodium nanoparticle of fig. 1, one obtains for n = 
1 (n = 2) the temperature T M /T F w 0.03 (0.04). This 
corresponds to an electronic temperature of Tm = HOOK 
(1600K). 

The curves in fig. 2 can be related to the experimentally 
accessible decay rate 



At \ T 



(13) 



In the case of an exponential relaxation, T t is indepen- 
dent of the time t. Since the differential transmission 
depends on time only implicitly through its dependence 
on temperature, the relaxation rate T t can be expressed 
as T t = -7T T d(lnT)/dt. The last factor on the right- 
hand side corresponds to the relaxation rate of the elec- 
tronic temperature and requires the precise knowledge of 
the time dependence T{t) of the temperature. While it 
determines the order of magnitude of T t , the dynamics 
of the relaxation process and in particular its frequency 
dependence is characterized by 
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This rate is independent of the time dependence T(t) of 
the temperature. 

In the sequel of the paper, we will use Tt taken at 
T = Tm in order to describe the relaxation process. This 
quantity is related to the slope of the differential transmis- 
sion at Tm, i.e., at the right end of the curves shown in 
fig. 2. The inverse of the relaxation parameter (14) is pre- 
sented in the inset of fig. 2 as a function of the frequency. 

For frequencies far away from the zeros of AT/T, the 
relaxation parameter is almost frequency independent and 
2(TmTt) 1 ~ 1. The factor of two arises since the differ- 
ential transmission is proportional to T 2 in this regime. In 
contrast, a strong frequency dependence appears close to 
lo\. For frequencies slightly below oj\ the initial slope of the 
differential transmission is reduced and eventually even 
changes sign. Therefore the inverse relaxation parame- 
ter exhibits a divergence. The increase of the differential 
transmission with decreasing temperature for the upper- 
most dashed curve (u/Up = 0.9849) can be traced back 
to the shift of uj\ which overcompensates the amplitude 
relaxation for this value of u> (see the inset of fig. 1). For 
frequencies slightly above uu\ we obtain very large slopes 
that decrease as we go further away from this frequency. 
As shown in the inset of fig. 2, between u\ and w\ a region 
appears where the relaxation parameter takes on values 
comparable to those far away from the resonance. While 
the anomaly around lo\ resembles the one close to w\, its 
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Fig. 2: The differential transmission normalised to the value 
at T/Tf = 0.05 is shown as a function of temperature. The 
system parameters are the same as in fig. 1. The solid lines 
correspond to the frequencies oj/^ip' = 0.9, 0.983 and 0.9847 
from the lower to the upper curve while the dashed lines cor- 
respond to u>/Wp> = 0.9852, 0.9854 and 0.9849 from the lower 
to the upper curve. In the inset, the inverse of the relaxation 
parameter (14) is shown as a function of the frequency for 
Tm/Tf — 0.05. The frequencies of the solid lines in the main 
graph lie on the left branch in the inset while the dashed lines 
lie on the middle branch. 

quantitative description relies on less justified hypothe- 
ses like the absence of an intcrband transition close to u)^ 
and the validity of the Breit-Wigner form over a broad 
frequency range. 

In the experiment of ref. [2] on copper nanoparti- 
cles, a differential transmission lifetime was observed at 
w/u4? ~ 0.99 which is almost a factor of two larger than 
the lifetimes found at lu/lo^ 0.97 and 1.025. Using 
the initial electronic temperature estimated to be about 
800 K [2] and the parameters of the experiment k F a = 68 
and ep/hujsp — 3.2, we find that the frequency at which 
the lifetime is a factor of two larger than far from the 
anomaly reproduces quantitatively the experimentally ob- 
served value. In addition, the increase of the lifetime is 
present only inside a very narrow frequency range whose 
width is smaller than the separations of the experimental 
points. Thus, in agreement with the experimental find- 
ings, a variation of the frequency by only 2% can yield an 
enlargement of the differential transmission lifetime by a 
factor of two. 

In fig. 3 we present T^ 1 for the sodium nanoparticle of 
figs. 1 and 2 as a function of the initial temperature Tvi 
for various frequencies around w^(Tm). Consistently with 
our findings for the frequency dependence of Tt, we find 
for frequencies far from o>£ (cf. the two lower curves) an 
almost frequency-independent relaxation parameter which 
only weakly increases with the maximum temperature Tm. 
On the other hand, for frequencies just below uj\ (cf. the 
uppermost curve), we see that T^ 1 is rapidly increasing 
with Tm beyond Tvi/Tp = 0.03. For higher initial tem- 
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Fig. 3: The inverse relaxation parameter (14) for the same 
system parameters as in fig. 1 is shown as a function of the 
maximum electronic temperature Tm in the nanoparticle. The 
frequencies are oj/J^ = 0.9, 0.983, 0.9847, and 0.9849 from 
the lower to the upper curve. 

peratures Tm, T^ 1 will eventually become negative as dis- 
cussed above. 

The results presented in fig. 3 can be compared with 
experiments because the initial temperature Tm is re- 
lated to the experimentally controlled pumping energy 
density which determines the energy initially deposited 
in a nanoparticle (sec eq. (12)). By measuring the relax- 
ation rate on and off resonance, it was found in ref. [4] 
that at an estimated average flux of one pump photon per 
nanoparticle a transition from a frequency-independent 
to a frequency-dependent relaxation rate occurs. Above 
this threshold, the inverse relaxation rate on resonance 
is significantly larger than off resonance. This scenario 
is clearly reproduced by the theoretical results shown in 
fig. 3. Consistently with the experimental findings, the 
lifetime can be increased by a factor of two when the 
number of absorbed photons is larger than one, i.e., for 
Tm/Tf > 0.03. 

While the slowdown of the transmission dynamics has 
been measured, an acceleration has not yet been reported 
in the literature. However, extracting the frequency- 
dependent relaxation times close to the resonance is quite 
difficult because there the initial value of AT/T used 
to normalise the data is very small. Systematic exper- 
iments measuring the differential transmission relaxation 
as a function of frequency would nevertheless be highly de- 
sirable in order to confirm unambiguously that the mech- 
anism we propose is at work. 

In summary, we have shown that the relaxation dynam- 
ics of the differential transmission of a nanoparticle ex- 
hibits an anomaly close to the surface plasmon resonance. 
This anomaly is due to the temperature dependencies of 
the position and the width of the surface plasmon reso- 
nance and manifests itself as a slowdown for frequencies 
below the plasmon frequency. In addition, in a narrow 
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frequency window we find a faster relaxation. 

The slowdown of the relaxation of the differential trans- 
mission close to the plasmon resonance has been observed 
experimentally for copper and silver nanoparticles [2,4,5], 
and several of the key findings are reproduced by our ap- 
proach, including the fact that the anomaly disappears 
under weak excitations or for relatively large clusters. 

The previously proposed explanation of the anomaly as 
a size-dependent many-body effect based on the resonant 
scattering of the d holes into the conduction band [7] can- 
not be invoked in the absence of interband transitions close 
to the surface plasmon resonance, as it is the case for silver 
nanoparticles. In contrast, the mechanism presented here 
does not rely on interband transitions and applies to a 
large variety of metallic nanoparticles. This demonstrates 
that the interband transition is not essential for the ap- 
pearance of the relaxation anomaly. In particular, it leads 
us to predict that a considerable slowdown of the trans- 
mission relaxation should also be observable in alkaline 
nanoparticles. 

Our explanation does not invoke any special behaviour 
of the electron dynamics at the plasmon frequency. In our 
view, the relaxation of the electronic temperature fully de- 
scribes the cooling of the electron gas on timescales which 
are larger than the lifetime of the surface plasmon. The 
time evolution of this temperature may depend on physi- 
cal parameters, like the size of the nanoparticle [6, 19] or 
the medium in which it is embedded [20] , but the decreas- 
ing electronic temperature is not expected to be influenced 
by the frequency at which the optical transmission is de- 
tected. The anomaly with a strong frequency dependence 
of the relaxation rate appears only when the differential 
optical transmission is considered. 

* * * 

We thank J.-Y. Bigot, V. Halte and P.-A. Hcrvicux for 
stimulating discussions. RAJ benefited from the hospi- 
tality of the Institut fur Physik der Universitat Augsburg 
during the final phase of this work. Financial support 
through the French-German PAI program Procope and 
from the European Union through the MCRTN program 
is gratefully acknowledged. 

REFERENCES 

[1] Tokizaki T. et al, Appl. Phys. Lett, 65 (1994) 941. 
[2] Bigot J.-Y., Merle J.-C, Cregut O. and Daunois 

A., Phys. Rev. Lett, 75 (1995) 4702. 
[3] Perner M. et al, Phys. Rev. Lett., 78 (1997) 2192. 
[4] Halte V. et al, Phys. Rev. B, 60 (1999) 11738. 
[5] Bigot J.-Y., Halte V., Merle J.-C. and Daunois A., 

Chem. Phys., 251 (2000) 181. 
[6] Del Fatti N. et al, Chem. Phys., 251 (2000) 215; Voisin 

C, Del Fatti N., Christofilos D. and Vallee ¥., J. 

Phys. Chem. B, 105 (2001) 2264. 



[7] Shahbazyan T. V., Perakis I. E. and Bigot J.-Y., 

Phys. Rev. Lett, 81 (1998) 3120; Shahbazyan T. V. 

and Perakis I. E., Phys. Rev. B, 60 (1999) 9090. 
[8] Weick G., Ingold G.-L., Jalabert R. A. and Wein- 
mann D., Phys. Rev. B, 74 (2006) 165421. 
[9] Hervieux P. -A., Benabbas A., Halte V. and Bigot 

J.-Y., Eur. Phys. J. D, 24 (2003) 185. 
[10] Eesley G. L., Phys. Rev. B, 33 (1986) 2144. 
[11] Weick G., Ph.D. Thesis, Universite Louis Pasteur 

Strasbourg and Universitat Augsburg (2006), URL: 

http://tel.ccsd.cnrs.fr/tel-00103438. 
[12] Gerchikov L. G., Guet C. and Ipatov A. N., Phys. 

Rev. A, 66 (2002) 053202. 
[13] de Heer W. A., Rev. Mod. Phys., 65 (1993) 611; Brack 

M., ibid., 65 (1993) 677. 
[14] Barma M. and Subrahmanyam V., J. Phys.: Condens. 

Matter, 1 (1989) 7681. 
[15] Yannouleas C. and Broglia R. A., Ann. Phys. (N. Y.), 

217 (1992) 105. 
[16] Molina R. A., Weinmann D. and Jalabert R. A., 

Phys. Rev. B, 65 (2002) 155427. 
[17] Weick G., Molina R. A., Weinmann D. and Jalabert 

R. A., Phys. Rev. B, 72 (2005) 115410. 
[18] Hervieux P. -A. and Bigot J.-Y., Phys. Rev. Lett, 92 

(2004) 197402. 

[19] Darugar Q., Qian W., El-Sayed M. A. and Pileni 
M.-P., J. Phys. Chem. B, 110 (2006) 143. 

[20] Link S. and El-Sayed M. A., J. Phys. Chem. B, 103 
(1999) 8410. 



p-6 



